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Abstract: We study the contributions to CP violation in right-handed sneutrino decays 
induced by soft supersymmetry-breaking gaugino masses including flavour effects and pay- 
ing special attention to the role of thermal corrections. Using a field-theoretical as well as 
a quantum mechanical approach we compute the CP asymmetries and we conclude that 
for all the soft-supersymmetry breaking sources of CP violation considered, an exact can- 
cellation between the asymmetries produced in the fermionic and bosonic channels occurs 
at T = up to second order in soft supersymmetry-breaking parameters. Once thermal 
effects are included the new sources of CP violation induced by supersymmetry-breaking 
gaugino masses can be sizeable and they can produce the observed baryon asymmetry for 
conventional values of the B parameter. 
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1. Introduction 

The discovery of neutrino oscillations makes Icptogcncsis a very attractive solution to the 
baryon asymmetry problem [1,2]. In the standard type I seesaw framework [3], the singlet 
heavy neutrinos liave lepton number violating Majorana masses and when decay out of 
equilibrium produce dynamically a lepton asymmetry which is partially converted into a 
baryon asymmetry due to fast sphaleron processes. 

For a hierarchical spectrum of right-handed neutrinos, successful leptogenesis requires 
generically quite heavy singlet neutrino masses [4], of order M > 2.4(0.4) x 10^ GeV for 
vanishing (thermal) initial neutrino densities [4,5] (although flavour effects [6-9] and/or 
extended scenarios [10,11] may affect this limit). Low-energy supersymmetry can be in- 
voked to naturally stabilize the hierarchy between this new scale and the electroweak one. 
This, however, introduces a certain conflict between the gravitino bound on the reheat 
temperature and the thermal production of right-handed neutrinos [12]. A way out of 
this conflict is provided by resonant leptogenesis [13]. In this scenario right-handed neu- 
trinos are nearly degenerate in mass which makes the self energy contributions to the CP 
asymmetries resonantly enhanced and allowing leptogenesis to be possible at much lower 
temperatures. 

Once supersymmetry has been introduced, leptogenesis is induced also in singlet sneu- 
trino decays. If supersymmetry is not broken, the order of magnitude of the asymmetry and 
the basic mechanism are the same as in the non-supersymmetric case. However, as shown 
in Refs. [14, 15], super symmetry-breaking terms can induce effects which are essentially 
different from the neutrino ones. In brief, soft supersymmetry-breaking terms involving 
the singlet sneutrinos remove the mass degeneracy between the two real sneutrino states of 
a single neutrino generation, and provide new sources of lepton number and CP violation. 
In this case, as for the case of resonant leptogenesis, it is the sneutrino self-energy con- 
tributions to the CP asymmetries which arc resonantly enhanced. As a consequence, the 
mixing between the sneutrino states can generate a sizable CP asymmetry in their decays. 
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This scenario was termed "soft leptogenesis" . Altogether it was found that the asymmetry 
is large for a right-handed neutrino mass scale relatively low, in the range 10^ — 10^ GeV, 
well below the reheat temperature limits, what solves the cosmological gravitino problem. 
However in order to generate enough asymmetry the lepton- violating soft bilinear coupling, 
B, responsible for the sneutrino mass splitting, has to be unconventionally small [14-17] *. 

In soft leptogenesis induced by CP violation in mixing as discussed above an exact 
cancellation occurs between the asymmetry produced in the fermionic and bosonic channels 
at T = 0. Thermal effects, thus, play a fundamental role in this mechanism: final-state 
Fermi blocking and Bose stimulation as well as effective masses for the particle excitations 
in the plasma break supersymmetry and effectively remove this degeneracy. 

In Rcf. [20] the possibility of soft leptogenesis generated by CP violation in right- 
handed sneutrino decay and in the interference of mixing and decay was considered. These 
new sources of CP violation (the so called "new ways to soft leptogenesis") are induced 
by vertex corrections due to gaugino soft supersymmetry-breaking masses. Some of these 
contributions, although suppressed by a loop factor and higher order in the supersymmetry- 
breaking parameters are relevant because they can be sizeable for natural values of the 
B parameter. Furthermore it was found that, unlike for CP violation in mixing, these 
contributions did not require thermal effects as they did not vanish at T = 0. 

In this work we revisit the role of thermal effects in soft leptogenesis due to CP vio- 
lation in right-handed sneutrino decays induced by gaugino soft supersymmetry-breaking 
masses. In Sec.2 we describe the one-generation see-saw model in the presence of the 
soft supersymmetry-breaking terms and compute the relevant CP asymmetries in a field- 
theoretical approach. We find that for all soft supersymmetry-breaking sources of CP 
violation considered, at T = the exact cancellation between the asymmetries produced 
in the fermionic and bosonic channels holds up to second order in soft supersymmetry- 
breaking parameters. In Sec. 3 we recompute the asymmetries using a quantum mechanical 
approach, based on an effective (non hermitic) Hamiltonian. We find the same T de- 
pendence of the resulting CP asymmetries. Finally in Sec. 4 we present our quantitative 
results and determine the region of parameters in which successful leptogenesis induced 
by the different contributions to the CP asymmetry is possible including the dominant 
thermal corrections as well as flavour- dependent effects associated with the charged lepton 
Yukawa couplings in this scenario. 

2. The CP Asymmetry: Field Theoretical Approach 

The supersymmetric see-saw model could be described by the superpotential: 

W = ^MijNiNj + YijE^pNi^H^, (2.1) 

where i, j = 1, 2, 3 are flavour indices and A^j, Lj, H are the chiral superfields for the right- 
handed (RH) neutrinos, the left-handed (LH) lepton doublets and the Higgs doublets with 
ea/3 = -e/3a and ei2 = +1. 

'Extended scenarios [18,19] may alleviate the unconventionally-small-B problem. 
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The relevant soft breaking terms involving the RH sneutrinos and SU(2) gauginos 
A2 are given by ^ 

jCsoft = - (^AijYije^pNii'^h'^ + ^BijMijNiNj + ^tti^XPl^ + h.c.^ . (2.2) 

The Lagrangian for interaction terms involving RH sneutrinos TVj, the RH neutrinos 
Ni and the A2 with (s)Ieptons and higgs(inos) can be written as: 

T PL{a^)ai\th^* - ^rPL(a3)a/3A0/i^*) + h.c. (2.3) 

ij = (fi,£^), = {viii^^ are the lepton and slepton doublets and, = {h^^hP) and 

= , h'^^ , are the Higgs and higgsino doublets. denotes A^ for aP = 01 and A^ 
for a(3 = 10 with cji^a being the Pauli matrices, and Pl,r are the left or right projection 
operator. 

The sneutrino and antisneutrino states mix with mass eigenvectors 



v2 

N_i = ^(e^*/2iVi - e-**/2^*), 
v2 



(2.4) 



where $ = arg(i?M) and with mass eigenvalues 



(2.5) 



Prom (2.3) and (2.4), we can write down the Lagrangian in the mass basis as 



1 

71 



A^4 



'^13 



h Pl^^ + {Aij + Mi)i'^hl^ 



\/2 



h PlP + {Aij - MiYqh^ 



h PLNit + NiPLifh^ 



-92 (xtPL{<7l)api?it - -^XPMaP^fit 

+T PL{oi)aP^h'^* - ^rPL{a3)^iXlhf^*^ +h.C.. 



(2.6) 



In what follows, we will consider a single generation of A'^ and N which we label as 1. 
We also assume proportionality of soft trilinear terms and drop the flavour indices for the 
coefficients A and B. 



^The effect of U{1) gauginos can be included in similar form. 
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As discussed in Refs. [14,15,20], in this case, after superficld rotations the Lagrangians 
(2.1) and (2.2) have two independent physical CP violating phases: 



which we choose to assign to A and to the gaugino coupling operators (the last two lines 
which are multiplied by g2 in Eq.(2.6) ) respectively. So for the calculations below we will 
take M, B, m2 and Yife to be positive real and A with phase 0a and define a complex 
coupling g2 = g2Q^p{i4>g) respectively. 

As discussed in Rcf. [15], when T ^ AM± = — M_, the two singlet sneutrino 
states are not well-scparatcd particles. In this case, the result for the asymmetry depends 
on how the initial state is prepared. In what follows we will assume that the sneutrinos are 
in a thermal bath with a thermalization time shorter than the typical oscillation times, 
AM^^, therefore coherence is lost and it is appropriate to compute the CP asymmetry in 
terms of the mass eigenstates Eq.(2.4). 

We compute the relevant decay amplitudes following the effective field-theoretical ap- 
proach described in [21], which takes into account the CP violation due to mixing and 
decay (as well as their interference) of nearly degenerate states by using resummed prop- 
agators for unstable mass eigenstate particles. The decay amplitude .A"* of the unstable 
external state iVj defined in Eq. (2.4) into a final state (a^ = s^, fk with = ^^/i^ and 
fk = ^k^^) described by a superposition of amplitudes with stable final states: 



cl>A = avg{AB*), 
^9 = -j^'^ziBml), 



(2.7) 
(2.8) 




(^5? +zV±'=^^^(M^)) - (a^'^ +iV^'=^'^^(M^)) 



,;vabs 



(2.9) 



M2 -M| + iSabs' 



.■^abs 





^ab are the absorptive parts of the iV{, — Na self-energies (see Fig. I): 




with 




M(FFt) 
Ait 



V = 



k 



11 



(2.15) 
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Figure I: Feynman diagrams contributing to the RH sneutrino self-energies at one-loop. 




Figure II: Feynman diagrams contributing to the RH sneutrino decay vertex at one-loop. 



yOfcabs absorptive parts of the vertex corrections (see Fig. 2): 



Stabs / 2\ 



ilfc 3772-2 2i "^2 

-(52) In 2_L,^2 ^«/^ 



V2 327r 

Yik 3m2 2 



-|- 777-2 



abs 



abs 



y) 



327r 

Yik 31712 



(52)' In 



_m2_ 

-|- 1712 



V2 327rp2 

. Ylk 3777-2 

' ^2 327rp' 



(^*+M)(5^)^ln 



p2 -|- 7772 



■{A* -M){~g*r In 



777-9 



p2 -|- 777-2 



[u{pe)PRv{ph)]eai3, 

[u{pf)PRv{ph)]eaf3 ■ 



(2.16) 
(2.17) 
(2.18) 
(2.19) 



2.1 The CP asymmetry 

The CP asymmetry produced in the decay of the states A'^j=-|- which enters into the Boltz- 



■''We notice that there is an irrelevant global i factor in the tree level and one- loop y^*"^^^ amplitudes 

compared to A"*" and \!°^^^^ arising from the particular choice of global phase in the definition of A''- in 
Eq.(2.4). 



- 5 - 



mann Equations (BE) is given by: 



i=±,ak 



(2.20) 



where we denote by 7 the thermal averaged rates. In the rest frame of N±, (2.20) simphfies 
to: 



-+k 



/ 

+k 



J 



where 



Sfc|2 



^±k 



(2.21) 



(2.22) 



i^^'^n 4VM± 



fk\2 \ Jk 



(2.23) 



In Eqs.(2.22) and (2.23) cl'^,c{'' are the phase-space factors of the scalar and fermionic 
channels, respectively. As long as we neglect the zero temperature lepton and slepton 
masses and small Yukawa couplings, these phase-space factors are flavour independent and 
they are the same for i = it. After including finite temperature effects in the approximation 
of decay at rest of the N± they are given by: 

c{{T) = cL{T) ^ cf{T) = (1 - X, - xj^)X{l,xe,xj^) [l - f^] \l - /f 1 , (2.24) 



4(r) = ci(r) = c^T) = X{l,xn,x^) [1 + f, 



h J 



eg 



where 



eg _ 



M exp[R./T]-l' 



■ ■ 

1 



(2.25) 

(2.26) 
(2.27) 



'M exp[£;^yr] + i' 

are the Boltzmann-Einstein and Fermi-Dirac equilibrium distributions, respectively, and 



\{l,x,y) = ^{l + x-yf-Ax, x„ = ^^. 

The thermal masses for the relevant supersymmetric degrees of freedom are [22]: 

'3 



ml{T) = 2m|(T) 
m}{T) = 2mjiT) 



(2.28) 
(2.29) 

(2.30) 
(2.31) 
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Here 52 and gy are gauge couplings and At is the top Yukawa, renormalized at the appro- 
priate high-energy scale. 

Substituting (2.9) and (2.10) into (2.22) and (2.23) we get in the numerators: 



-4 I -Im 
-l-Im 



^afc*^afc^abs 



(M^-M|)2 + |S|bs|2 



^afc*yafcabs^^2^ 



-|-Im 



ya^abs* ^jy^2^-j^afc j^abs _ ^cifc*ya.fcabs^^2^-jj^abs 



-labs 



(2.32) 



(M^-M|)2 + |sabs|2 J ' 



where we have used the relations T,^^ = T,^^ and S^^^* 
terms of order 6g and higher are ignored with 

^ ~ M' M' M' 



■^abs 



The ~ sign means that 



(2.33) 



The three lines in (2.32) correspond respectively to (i) CP violation in N mixing from the 
off-diagonal one-loop self-energies (this corresponds to the effects originally considered in 
Refs. [14, 15]), (ii) CP violation due to the gaugino- mediated one-loop vertex corrections 
to the N decay, and (iii) CP violation in the interference of vertex and self-energies. 

In the denominator of (2.22) and (2.23) we consider only the tree-level amplitudes 
+ = 2\Al^\\ with = Y^l [|^|2 ± 2MRe(yl)] and \A^^\^ = Y-^^^l. 

Using the explicit forms in Eqs.(2.14) and (2.19) we find that up to order 5l, the three 
contributions to the CP asymmetry from scalar and fermion decays verify: 

c/(T) 



.5 _ c-(T) s 



sV 
^±k 



^±k 



c'{T)+cf{T) 
(f{T)+cf{T) 
c'{T) + cf{T) 



^±k — 



^±k — 



^±k 



c'{T) + cf{T) 
c'{T) + cf{T) 
c'{T) + cf{T) 



^±ki 



with 



^±k 



-Kl^ ( 1 T sin 



M 



2BT 



2M J 



4^2 + r2 ' 



SK2 



"8 H 



In- 



ml 



mla2 m2 \A\ 



mi 



In- 



-hM2 



mt, 



4 MM mi + M"^ 



\A\ B 

— sin {(f)A + 2(t>g) - — sin (2(^g) ± sin (2(^g) 



sin ((^a) cos {24) g) 



^' 452 + r2 ' 



(2.34) 

(2.35) 
j(2.36) 
(2.37) 



2 

where we have defined a2 = || and the flavour projections Ky. 



\Yik? 



(2.38) 
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Summing up the contribution from the decays of and one gets the three 
contributions to the CP asymmetry in Eq. (2.20) 

= KlABF{T)f, (2.39) 



i^sin(0A + 2.^g)--^sin(205) 



Abf{T), 



4 M ml + W^ 
= KlABF{T)e'' , (2.40) 

4(T) = ^i^ln-^sm(,.)cos(2,,)^A.HT), 

= KI/^bf{T) e', (2.41) 

where 

The asymmetry in Eq.(2.39) is the contribution to the lepton asymmetry due to CP 
violation in RH sneutrino mixing discussed in Refs. [14-16]. Eqs.(2.40) and (2.41) give 
the contribution to the lepton asymmetry related to CP violation in decay and in the 
interference of mixing and decay. They have similar parametric dependences as the ones 
derived in Ref. [20]. However as explicitly shown in the equations, the scalar and fermionic 
CP asymmetries, (2.34), cancel each other at zero temperature. Consequently we find that, 
up to second order in the soft supersymmetry-breaking parameters, all contributions to the 
lepton asymmetry in the soft supersymmetry scenario require thermal effects in order to 
be significant. 

We finish by noticing that in this derivation we have neglected thermal corrections 
to the CP asymmetry from the loops, i.e., we have computed the imaginary part of the 
one-loop graphs using Cutkosky's cutting rules at T = 0. 

3. CP Asymmetry in Quantum Mechanics 

In this section wc recompute the asymmetry using a quantum mechanical (QM) approach, 
based on an effective (non hermitic) Hamiltonian [14,15,20]. In this language an analogy 
can be drawn between the N-W system and the system of neutral mesons such as K^-K 
and its time evolution is determined in the non-relativistic limit by the Hamiltonian: 




H= \ B L\--Ata h (3-1) 



with r given in Eq.(2.15). 

In Refs. [14, 15,20] the QM formalism was applied for weak initial states N and W . 
In practice it is possible to use the formalism to study the evolution of either initially 
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weak or mass eigenstates. So in order to study the dependence of the results on the choice 
of physical initial conditions we will compute the asymmetry in this formalism assuming 
either of the two possibilities for initial states. So we define the basis: 



iVi 



7V2 = e'l^ (bN - aiVt) . 



(3.2) 



The mass basis, Eq.(2.4) corresponds to {a,b,P) = (;^!;^'~f)- Assuming that the 
physical initial states were pure N and corresponds to {a,h,(3) = (1,0, tt). 

The decay amplitudes of Ni and N2 into fermions fk = (■%h'^ including the one-loop 
contribution from gaugino exchange are: 



4' 



Ai 



Yika - tttS (6M + aA) (52)' ^If [u {Ph) Plv {pi)\e^, 



2M2 



167r 



(6M - aA*) {glf ^/^| [u{p,) Prv (pn] e^, 



31u 
2M2 



(3.3) 



where the A denotes the decay amplitudes into antifermions. The corresponding decay 
amplitudes into scalar = i%h'^ and 



where 
Re(J/) = fR = 

Re(J5) = SR = 



A\^ 
If 



Al' 



Al' 



1 

vr 

1 

TT 



Fife (aM + hA*) - {hf \ eed, 

(6M + aA) - ^a {'9*2)' \ ^cd. 



yifc (6M - aA*) + (52)' ^4 [> fed, 

Fife (6A - aM) - ^5 (5^)2 [> 



- ( In- 



+ M2 



+ Li2 



+ M2 



In- 



+ M2 



+ Li2 



mi 



Im(7/) = ji = Im{Is) = SI 



-In- 



rrin 



771-2 + ' 

The eigenvectors of the Hamiltonian in terms of the states A^^i and N2 are: 

^2 



(3.4) 



-C(2) 

-C(2) + Bo(M2,m2,0) +So(M2,0,m2) 

(3.5) 



Nl) = {ap + bq) Ni)+ e-'l^ {bp - aq) 
Nh) = (ap - bq) Ni)+ e''^ {bp + aq) 



Nn 



(3.6) 
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where 



q 

P 



-1 -341 sin (<^^) 



■cos {(pA) - -- 



BM—y-r^j ^^^^ 2M252 

At the time t the states Ni and iV2 have evolved into 

Ni,2{t)) = ^ {[eUt) + eH{t) ± Co {eL{t) - enim 



sin {2(f) a) 



Nl,2 



(3.7) 



where 



Co = a6 ( ^ + ^ 
q p 



Ci = a -, C2 = a -, 



P 



P 



and 



The total time integrated lepton asymmetry is 

QM _ t=l,2,afc 

i=l,2,ak,k 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



where T{Ni at) are the time integrated decay rates which from Eq.(3.8) are found to be 



TiNi ^ ak) 



1 c"-k 



4167rM 



Gi+ + 



+2 



Re ( ) G,? - Im ( A^^A^i, ) Ci 



. (3.12) 



T{Ni (Ik) can be obtained from Eq.(3.12) with the replacement ^d"* — >■ A'^*' . We have 
defined the time integrated projections 



Gl,2+ = 



2 

±8 



1 1 

+ 



1 — y I + 



+2|Cor 



1 



1 



y 



Re (Co) -^-Im (Co) ^ 
I — y 1 + 



1 — 1 + 



(3.13) 



'^11(22) 



^11(22) 



2|Cl,2 

2<fRc 



1 - y2 i + x' 



(i^'^C- 



1(2) 



l-y2 



Im 



3^*'^C 



1(2) 



1 



±2Re 
2 ilm 



1 - y2 1 + 



1(2) 



1 - r 



+ Re 



±2Im 



eT^^Co*Ci(2) 



1 



1 +X2 



1 — 1 + x^ 



(3.14) 



(3.15) 
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in terms of masses and width differences coefficients 



X = 



V 



B 
T 



IT\A[' 



2 BM2 



sin^ {(I)a) , 



2T 



= — cos {(Pa) • 



Substituting Eqs.(3.12)-(3.15) one can write the numerator in Eq.(3.11) as 



(3.16) 
(3.17) 

(3.18) 



with 



1 c"'= 



,r- + y- 



2 167rM(l-y2)(i + a;2) 



{n\H\AT\ 



^2 



(|Ci|'-|C2 



akl"^ 



(3.19) 



+2 
-2 



Re (^Al>'*Al'' - If^^^'=) Re (e-'f^C^d^ - Re + Re (e^^CfJCa' 

Im - ]4f^*]4f ) Im (^e'^'^Co^Ci) - Im + If *]4f^) Im (^e^^C^Ca)] }, 



r'^fc 1 r / _ 2 _ 2\ 



+ ( l^f 1^ 



+ l^2l 



^2 



(3.20) 



+y 
-y 



Re (y4^'=*A^'= - Af*Af^ Re (e'^^Ci) + Re [aI^^AI'^ - AI>'* Af^ Re (e'^C2)] 
Im ( A\^*Al^ - If^*!^) Im [e-'<^C^ + Im - If *Af^) Im (e^'3C2)] }, 



(3.21) 



c"-^ r ( / o I — ^ 2 „ — ^ 2\ 

= 16^M(1 + X2) { " ^^""^^"^ ('"^l"' " " + ^2*= J 

Re - If *If ) Re (e'^'^Ci) + Re (aS^^* A"^ - If *If ) Re (e'<^C^ 

Im (a^**^^''- - If *If ) Im (e-^^Ci) + Im (aI^"" A\^ - If *If ) Im (e}f^C^ } . 

In writing the above equations we have classified the contributions as resonant, (non- 
resonant) , R (NR.) , depending on whether they present an overall factor ^ipp- (or no 
at all). We have labeled the remainder as interference term, /. 

After substituting the explicit values for the amplitudes and the coefficients and ne- 
glecting all those terms which cancel in both basis we get that 



R 

k ) 



Ar/fc.^« = -c^Arf 



^Y^k,NR = c'^Ar^^, 



(3.22) 



^We use the expression of Th —Yl from Ref. [20]. Notice that with this definition Th —^l 7^r+ — r_ 
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with 



1 

in 
3 
167r 



Y^^a2 In 



6^)2 + (2a6)2 cos 2/3] |^| sin((/.A)-7^ 



m2 



1 



+yM (2(2a6)^ + (a 



+ M2 M 1 - 7/2 

2 



''2 

\2 



6^)2 COSI 



167r ^'^"^ mi + M2 M 1 + x2 



(l-2/2)(l + ^2) 

[-\A\sm{(f)A + 2(t)g) 

2(3)) sm(203)] , 

1^1 sm.{<j)A) cos(2/3) cos(2^g). 



(3.23) 

(3.24) 
(3.25) 



Eq.(3.22) explicitly displays the cancellation of the asymmetries at T = also in this 
formalism for either initial mass or weak eigenstate right-handed sneutrinos. ^ 

Introducing the explicit values for the coefficients for initial weak RH sneutrinos and 
the expressions for x and y and expanding at order (5| we get 



NR,QMw 



(T) 



M 

3K9a2 1712 



In 



S2 + r2 



4 M ml + M'^ 



\A\ 

— sin(</)A)cos(2(/)g) 



(3.26) 



+ 



B 



sin (2^g 



4'«^^(r) = 



3K^^a2 7112 \A\ 



In- 



ml 



4 MM m'i + M'^ 



2M 

sin ((?!)a) cos {2(j)g) 



r2 



fi2 + r2 



Correspondingly for initial N± states one gets 



R,QMm 



(T) = K[ 



\A\ 



BY 



sm I 



NR,QMm,rj,^ ^ 3.R:gQ!2 m2 , 

fe ^ > AM 



B2 + r2 

2 

m-2 



Asf(T), 



(3.27) 
(T) .(3.28) 

(3.29) 



m. 



+ M2 



1^1 



sin(^A) cos(2^g) 



B 



sin (20p) 



J,QM,n^^ ^ 3i^fcQi2 m2 1^1 



4 MM + M2 



2M 

sin (0^)cos(2^g) 



r2 
52 + r2 



(3.30) 

ABi.(r)(.3.31) 



Comparing Eqs.(3.29)-(3.31) with Eqs.(3.26)-(3.28) and Eqs.(2.39)-(2.41) we find that 
they show very similar parametric dependence though there are some differences in the nu- 
merical coefficients. In particular we find that e^'^^ , e^''^''^ and the i?-dependent (second 
term) in either the weak or mass basis ^^^^^^ coincide with e^, ej, and the B-depcndent 
term in derived in the previous section with the redefinition A 2A, B 2B and 
sm{(f)A) — ±sin(^^). We find only some differences in the phase combination which ap- 
pears in the B independent term in the asymmetries e^'*^^ and as seen in Eqs.(2.40) 
and Eqs.(3.27) and (3.30). 



^We have traced the discrepancy with Ref. [20] to a missing cos{(pf — 0s) factor in their expression for 

sin Ss in their Eq.(19). Once that factor is included, £2* in their Eq.(22) cancels against e^* in their Eq.(25), 
and e* in their Eq.(23) cancels against e** in their Eq.(24) so that the total asymmetry is zero at T = 0. 
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4. Results 



Next we quantify the parameters for which successful leptogenesis induced by the differ- 
ent sources of CP violation discussed in the previous sections is possible by solving the 
corresponding set of Boltzmann Equations (BE). 

The relevant classical BE describing the decay, inverse decay and scattering processes 
involving the sneutrino states in the framework of flavoured soft leptogenesis were derived 
in detail in Ref. [16] assuming that the physically relevant sneutrino states were the mass 
eigenstates (2.4). As we have seen in the previous section, the choice of the physical basis 
for the sneutrino states docs not lead to important differences in the parametric form of 
the generated asymmetries, thus in what follows we will present our results assuming that 
the relevant sneutrino states are the mass eigenstates. 

Including the N± and N decay and inverse decay processes as well as all the AL = 1 
scattering processes induced by the top Yukawa coupling the final set of BE takes the form 
II : 

sHz'-^ = - - l) (t. + + 47« + 47?) + 2,i^^ + 4,(^)) , (4.1) 

^) + if + 37.2 + 2,i'^ + 27f ) + 2jP + r + 27f )) 
+lNE'kiT)^, (4.2) 



^tot _ _ ( ^tot 



sHz — — ^ = — < 
dz 



N 




' 7I2 ) ■ (4.3) 



eg 



In the equations above Yc'^ = and Y'^^T > M) = 90C(3)/(47r^c/*), where 5* is the 

total number of entropic degrees of freedom, g* = 228.75 in the MSSM. 

Before proceeding let's point out that, as recently discussed in Refs. [23], for resonant 
scenarios, the use of quantum BE [23, 24] (QBE) may be relevant. For standard resonant 



"AL = 1 scattering involving gauge bosons, AL = 2 off-shell scattering processes involving the pole- 
subtracted s-channel and the u and t-channel, as well as the the L conserving processes from N and N pair 
creation and annihilation have not been included. The reaction rates for these processes are quartic in the 
Yukawa couplings and can be safely neglected. 
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leptogenesis they induce a T dependence in the CP asymmetry which can enhance the 
produced baryon number. However, as shown in Ref. [17] for soft leptogenesis from CP 
violation in mixing, due to the thermal nature of the mechanism already at the classical 
level, the introduction of quantum effects does not lead to such enhancement and the 
required values of the parameters to achieve successful leptogenesis are not substantially 
modified. Thus in this work we study the impact of CP violation in sneutrino decay and 
the interference of mixing and decay on the values of the parameters required for successful 
leptogenesis in the context of the classical BE as described above. 

In writing the BE relevant in the regime in which flavours have to be considered 
[6,7,11,25], it is most appropriate to follow the evolution of Y^^. where = ^ — YLk^ ~ 
Yl^^ = §-Y^k . This is so because A^. is conserved by sphalerons and by other MSSM 
interactions. In particular, notice that the MSSM processes enforce the equality of fermionic 
and scalar lepton asymmetries of the same flavour. 

In Eqs.(4.1)-(4.3) we have accounted for the CP asymmetries in the N± two body 
decays to the order described in the previous sections. We have neglected higher order terms 
in supersymmetry-breaking parameters which could lead to differences in the distribution 
and thermal widths of N± and correspondingly we have written a unique BE for = 
Y~ + Y~ 

In Eq. (4.3) we have defined the flavoured thermal widths 

7^ = ^.%^, l?' = KhP, (4.4) 

where the different 7's are the thermal widths for the following processes (in all cases a 
sum over the CP conjugate final states is implicit): 

7^ = 7^ + = 7(iV± ^ hi) + 7(iV± ^ hi), 

if = 7{N± ^ tuq) , 

722 = -f{N±i ^ uq) = -/{N±q* i*u) = -f(N±u* ^ i*q), 
7^^) = 7(iV±£ ^ qu) = -i{N^l ^ qu) , 
7^ ) = ^(iV±u ^ Iq) = -f{N±q* ^ lu) , 
7^ ) = -f{N±q ^ lu) = -i{N±u ^ Iq), 
7^ = l{N±t ^ qu), 

= -f{N±q ^ lu) = -f{N±u ^ iq), 
7^ = 7(iV ^ ih) + 7(iV ^ i*h), 
7^°) = 7(iV| ^ qu) = ^{Ni ^ qu), 
7^^) = j{Nq ^ tu) = 7(iV ^ tq) , 
7J2) = 7(ivr ^tq)= j{Nq* ^tu) , 
^(3) = ^(iV£ ^ qu) , 

7^ ) = 7(iV ^ Iq) = -f{Nq ^ lu) . (4.5) 

The explicit expressions for the 7's in Eq. (4.5) can be found, for example, in [27] for 
the case of Boltzmann-Maxwell distribution functions and neglecting Pauli-blocking and 
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stimulated emission as well as the relative motion of the particles with respect to the 
plasma. 

The value of Aap depends on which processes are in thermal equilibrium when leptoge- 
nesis is taking place. As we will see below for any of the considered sources of CP violation, 
the relevant temperature window around M ~ T corresponds to T < (1+tan^ /?) x lO^GeV. 
In this regime the processes mediated by all the three charged lepton (e, r) Yukawa cou- 
plings are in equilibrium i.e. they are faster than the processes involving N± and one 
gets [26] 

\ 

A=\ I -¥o ■ (4-6) 




After conversion by the sphaleron transitions, the final baryon asymmetry is given by 

k 

where '^^Y^i^^z — ^ oo) can be parametrized as 
k 

^ Ya, {z^oo) = -2(e-^ + + e') r?fla , (4.8) 
k 

with e'^, e^, and are given in Eqs.(2.39)-(2.41). rysa is the dilution factor which takes 
into account the possible inefficiency in the production of the singlet sneutrinos, the erasure 
of the generated asymmetry by L-violating scattering processes and the temperature and 
flavour dependence of the CP asymmetry. It is obtained by solving the array of BE above. 
Within our approximations for the thermal widths, 77fla depends on the flavour projections 
K^, on the Yukawa couplings {YY^)ii and on the heavy mass M, with the dominant 
dependence on these last two arising in the combination 

{YY^)nvl = meffM, (4.9) 

where Vu is the vacuum expectation value of the up-type Higgs doublet, Vu = v sin/3 (t;=174 
GeV) . There is a residual dependence on M due to the running of the top Yukawa coupling 
as well as the thermal effects included in Abf although it is very mild. Also, as long as 
tan P is not very close to one, the dominant dependence on tan P arises via Vu as given in 
Eq. (4.9) and it is therefore very mild. 

In Fig. Ill we plot |77fla| as a function of rueff for M = 10'' GeV and for the the equally 
distributed flavour composition = K2 = = 1/3. We consider two different initial 
conditions for the sneutrino abundance. In one case, one assumes that the N population 
is created by their Yukawa interactions with the thermal plasma, and set Yf^{z ^ 0) = 0. 
The other case corresponds to an initial abundance equal to the thermal one, Yjy{z 



N ^ 

depending on the value of rrieff while for thermal initial conditions, on the contrary, 77 > 0. 



0) = Y^{z ^0). As discussed in Ref. [16] for zero initial conditions, r] can take both signs 
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10 ^ 10 ^ 10 10 10 ^ 10 ^ 

nieff (eV) 

Figure III: Efficiency factor \r]{ia.\ as a function of We// for M = 10'' GeV and tan/3 = 30 and for 
Ki = K2 = K^^ — 1/3. Tire two curves correspond to vanishing initial A'' abundance (solid black 
curve) and thermal initial N abundance, (dashed red curve). 

Introducing the resulting r^fla in Eqs.(4.8) and (4.7) we can easily quantify the allowed 
ranges of parameters for which enough asymmetry, Yb > 8.54 x 10~^^ [28], is generated. 
We plot in Fig. IV the resulting ranges for B and rrieff for the equally distributed flavour 
composition = K2 = = 1/3 and for \A\ = m2 = 1 TeV and tan/3 = 30. For different 
values of the CP phases and M as explicitly given in the figure. 

The upper panels in Fig. IV give the parameters regions for which the CP violation 
from pure mixing effects, e"^, can produce the observed asymmetry as previously described 
in [14-16]. Due to the resonant nature of this contribution, these effects are only large 
enough for B ~ 0{T) which leads to the well-known condition of the unconventionally 
small values of B and to the upper bound M < 10^ GeV. 

The central panels of Fig. IV give the corresponding regions for which CP violation from 
gaugino-induced vertex effects, e^, can produce the observed baryon asymmetry. Despite 
being higher order in 63 and including a loop suppression factor, 02, this contribution can 
be relevant because it is dominant for conventional values of the B parameter. However, 
in order to overcome the loop and 6s suppressions this contribution can only be sizeable 
for lighter values of the RH sneutrino masses M < 10^ GeV (within the approximation 
used in this work: 63 ^ 1, |j4|,m2 ~ ©(TeV)). 

The parameters chosen in the figure are such that the second term in Eq.(2.40) dom- 
inates so that the allowed region depicts a lower bound on B. Conversely, when the first 
term in Eq.(2.40) dominates, becomes independent of B. In this case, for a given value 
of M and 5s the produced baryon asymmetry can be sizeable within the range of m^ff 



values for which r/ga is large enough. For example for M = 10^ GeV, and m2 = |A| = 1 
TeV and | sin(^^ + 2^^) = 1| with vanishing initial conditions 

eV eV 

where each range corresponds to a sign of the CP phase sin(0yi + 2(j)g) 

Finally we show in the lower panels of Fig. IV the values of B and TTT-e// for which 
enough baryon asymmetry can be generated from the interference of mixing and vertex 
corrections e^, Eq.(2.41). Generically , is subdominant to e*^ since both involve the same 
CP phase sin((^A) while has additional 5s and loop suppressions: 

—3 mo , mo , , , F , . . 

Consequently as seen in the above equation and illustrated in the figure, can only 
dominate for extremely low values oi B {B <^ V) for which it becomes independent of B. 
Also we notice that for M < 10^ GeV and mg// > 10~^ eV the resulting baryon asymmetry 
generated by this contribution becomes independent of me// since the 'rn'^j:j, dependence 
from F^ cancels the approximate l/m^jj dependence of r/fja in this strong washout regime. 

In summary in this work wc have quantified in detail the contributions to CP violation 
in right-handed sneutrino decays induced by soft supersymmetry-breaking gaugino masses 
paying special attention to the role of thermal effects. Using a field-theoretical as well as 
a quantum mechanical approach we conclude that for all the soft supersymmetry-breaking 
sources of CP violation considered, an exact cancellation between the asymmetries pro- 
duced in the fermionic and bosonic channels occurs at T = up to second order in soft 
supersymmetry-breaking parameters. However, once thermal effects are included the new 
sources of CP violation induced by soft supersymmetry-breaking gaugino masses can be 
sizeable and they can produce the observed baryon asymmetry for conventional values of 
the B parameter. 
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Figure IV: B,meff regions in wliicli successful soft leptogenesis can be achieved when flavour 
effects are included with Kf = K2 = = 1/3 and for different sources of CP violation. In all 
cases we take A\ = m2 = 10^ GeV and tan/3 = 30 and different values of M and (t>A and 4)g as 
labeled in the figure (see text for details). The left (right) panels correspond to vanishing (thermal) 
initial N abundance . 
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